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ABSTRACT 

The ionospheric differential refraction, at distances of tens 
of kilometers or more, is one of the main problems 
affecting the capability of instantaneous carrier phase 
ambiguity resolution. And, as a consequence, the 
ionosphere affects the feasibility of centimeter- accuracy 
navigation with dual -frequency Global Navigation 
Satellite Systems such as present-day GPS, as well as 
future systems with three frequencies, like GALILEO and 
the Modernized GPS. We have shown in previous work 
thai a real-time tomographic model of the ionosphere, 
computed from the GPS reference network data, is 
accurate enough to allow GPS real-time carrier phase 
ambiguity resolution for a rover at hundreds of kilometers 
from the nearest reference site. In this work, we study the 
extension of such techniques to the improvement of 
instantaneous three-carrier phase ambiguities resolution 



algorithms, such as TCAR, at medium and long distances 
(from lens to hundreds of kilometers), and with a 
minimum of geodetic computation. The data sets used 
have been generated with a modified GNSS satellite signal 
generator. In particular, improvements in the instantaneous 
ambiguity resolution success rate with three frequencies 
can be attained from about 60% of all attempts in previous 
work to about 100% with the proposed technique, at 
distances of about 60 km and with low ionospheric values, 
and to more than 90% at distances of more than 100 km 
under Solar Maximum conditions. 

INTRODUCTION 

Long-range, real-time, carrier-phase DGPS, with sub- 
decimeter precision and separations of up to several 
hundreds of kilometers between the rover and the nearest 
station, requires the estimation or resolution of the carrier- 
phase ambiguities and other unknowns to achieve high 
precision. Effects that tend to cancel out between receivers 
in short-baseline DGPS, remain quite significant in long- 
range DGPSs: Broadcast orbit errors, tropospheric 
refraction correction errors, and floated ionospheric- free 
combination (Lc) biases. These additional unknowns have 
to be estimated simultaneously with the position of the 
user's rover, in the navigation Kalman filter. This filter 
needs to receive and assimilate GPS data for a while, until 
it has enough to make precise estimates of all the 
unknowns. Until then, (he position cannot reach the 
desired accuracy. 'Hi is could take the better part of one 
hour, and is known as the convergence period of the filter. 
The main reason for this delay is the estimation of the 
lloated Le biases. If one could resolve the corresponding 
ambiguities and then use the result to eliminate the biases, 
the convergence period would be much shorter, and 
precise navigation could begin that much sooner (Colombo 
et al., 1999, 2000). This is the main reason for attempting 
to resolve ambiguities in long-range navigation. 

The planned three-carrier satellite navigation systems offer 
the potential advantages of high success and integrity in 
instantaneous ambiguity resolution, with a minimum 
geodetic computation load. But it can be seriously affected 
by the ionospheric refraction, as is detailed below. 
Therefore, it makes sense to introduce an accurate 
modelling of the ionosphere in real time, as we have 
demonstrated in ( Hernandez- Pajares et al. 1999a, 2000 a, 
Colombo etal. 1999,2000). 
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Figure 1: Global Toial Electron Content map in which the 
vertical ionospheric delay, computed from actual GPS data, is 
represented for the target day of the experiment analyzed in this 
work, 17 March 2002 (units: I0" 1 TECU = 1.6 cm in LI). 

Several techniques have been proposed so far to improve 
the instantaneous positioning by using three- frequency 
systems, such as Three-Carrier Ambiguity Resolution, or 
TCAR, and Cascade Integer Resolution, or CIR (see for 
instance Vollath et al. 1998, Jung et al. 2000). They share 
a similar basic approach: the double differenced integer 
ambiguities are successively solved from the longest to the 
shortest beat- wavelength, including "extra- wide lane" and 
"wide lane" combinations of carrier phases (with 
wavelengths of 7.480m and 0.862m, respectively, in the 
datasets used in this work), and the first carrier LI (with a 
wavelength of 0. 1 90m). 

But, in particular, TCAR performance is strongly affected 
by the ionospheric refraction decorrelalion with distance. 
Indeed, as we will see, the ionosphere delay is a problem 
when its double differenced value is greater than 0.26 
Total Electron Content Units (TECU) 1 (about 4 cm in LI). 
And typical values can exceed this threshold, as it is 
shown in the vertical ionospheric Total Electron Content 
(TEC) map computed from GPS data in figure I. As a 
consecuence, an Integrated TCAR approach (ITCAR) has 
been developed by several authors (Vollath et al. 1998), 
including search algorithms and a navigation filter in 
which the ambiguities are part of the outputs and the 
ionospheric residual errors are coarsely estimated. But the 
ITCAR algorithm is still affected by the hick of knowledge 
of the double difference ionospheric refraction, limiting 
successful ambiguity resolution to cases where the 
separation between roving and fixed receivers is limited to 
a few tens of kilometers or more (Vollath el al. 2001). 

In this paper, we extend our previous work on the 
application of real-time ionospheric corrections to 
ambiguity resolution. These corrections are generated in 
real lime with a Kalman filler lhat assimilates dual- 
frequency carrier phase data from the network of reference 
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stations. They are used to eliminate the effect of the 
ionosphere from the double-diferenced data, making it 
possible to resolve the ambiguities in real time, with 
dual -frequency systems such as GPS. One example is our 
Wide Area Real Time Kinematic (WARTK) method 
(Hernandez-Pajarcs ct al. 2000, 2001). These accurale 
ionospheric corrections can be broadcast to the users. In 
turn, these correct and then process their "ionosphere free" 
data to resolve their ambiguities wilh a simple algorithm 
such as TCAR . As we explain in this work, this approach 
(WARTK-3) produces a better performance at long 
distances from the nearest reference site than ITCAR, 
allowing the instantaneous resolution of more than 90% of 
ambiguities over baselines of more than 100 km. 
Moreover, WARTK-3 is a quite simple algorithm wilh a 
low computation load for ihe user, compared with ITCAR. 
Also WARTK-3 is simpler than the dual-frequency 
algorithm WARTK, lhal we shall explain briefly in the 
next section. 

WARTK: PREVIOUS WORK WITH GPS 

We have shown in previous work lhat a real-time 
lomographic model of the Ionosphere is accurale enough 
to solve On-The-Fly the GPS carrier phase ambiguities in 
a rover at hundreds of kilometers from the nearest 
reference site, using the corresponding reference network 
corrections (figure 2), generated with a combined 
ionospheric and geodetic filter (see Hernandez-Pajarcs et 
al. 2000, 2001 for details). We have tried this form of 
ambiguity resolution for precise navigation, with sub- 
decimeter precision, and for making instantaneous 
atmospheric water- vapour determinations for 

meteorological applications, within wide-area networks 
with stations hundreds of kilometers apart. This approach 
fills the gap in the typical accuracy versus baseline length 
plot (see figure 3), between the Local Area DGPS (or 
GBAS) systems, with centimeter accuracy, but limited to 
distances below 20 km, and the new Global DGPS systems 
with sub- meter accuracy over thousands of kilometers 
(Bar-Sever et al., 2001). 

Long-Baseline (hundreds Km) OTF Ambiguity Resolution 




Figure 2: The Wide Area Differential GPS (WARTK) 
corrections, in particular the ionospheric ones, should be 
broadcast to the users. 
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WARTK has been tested in different scenarios: under 
hard ionospheric conditions typical of Solar Maximum, 
during Ionospheric Storms, and in the presence of 
Travelling Ionospheric Disturbances /TTO's (Hernandez- 
Fajares el al., 2001b). On the other hand, in the case of 
distances between stations of more than several hundred 
kilometers, strong electron content gradients can limit the 
performance of the technique. For this reason, an 
extension of the algorithm has been developed recently to 
include the case of reference stations separated thousands 
of kilometers (Hern£ndez-Pajares et al. 2002). 

ORIGINAL TCAR ALGORITHM FOR 3 
FREQUENCY SYSTEMS 

The planned three- frequency systems should offer more 
simultaneous observations, improving the chances of 
instantaneous (single epoch) ambiguity resolution, 
compared with the present dual- frequency GPS system. 
How ihis may work can be understood by studying one 
of the first algorithms proposed for 3 frequency navigation 
systems, the Three Carrier Ambiguity Resolution, or 
TCAR (Harris 1997). It consists of three basic steps: 
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Figure 3: Accuracy versus baseline length (from the nearest 
reference site) plot for different GPS differential systems. The 
WARTK approach tills a gap, with subdecimeter errors at 
hundreds of kilometers baselines. 

Step 1. Solve the extra-wide lane ambiguity (with a 
wavelength of -7.45 meters in our dataset) by adding a 
pseudo-range combination. Although in some eases the 
pscudo- range multipath can diminish the chances for 
succes, this error is typically small compared with the long 
wavelength of the extra- wide lane, and can not present 
much of a problem. 

The llrst step of the TCAR approach (see for example 
Vollaih el al. 1998) is to estimate the double differenced 
(between pairs of receivers and satellites) 
ambiguities VAN ilw of the double differenced extra-wide 
lane carrier phases VAL <<U . , using a combination of 
pseudo-ranges (or codes) that share the same value 
and sign of ionospheric delay as the wide lanes. . This is 



possible due to the. very long wavelength of the extra- 
wide lane combination, such the 7.45m in the case of the 
datasets generated for this study (Laboratory exp. 2000, 
and AGGA validation reports, 2000). 

The corresponding equations are: 

J3 J I 

= P'+ *>Jb„ +«„,/ + "»„ + £«+■■ 0) 

L x being the carrier phase observation (in length units) 
of frequency f x and wavelenght X x (see table 1 below), 
and P x the corresponding pseudorange observation. The 
multipath and observation errors for carrier phases and 
pseudoranges are indicated by m x , M x and E x , E x 
respectively (see maximum values of multipath and typical 
measurement errors in table 1). The indifferenced carrier 
phase ambiguity & x b x , that contains instrumental delays 
and, after double-differencing, becomes an integer 
multiplier value of the wavelength, X X VAN X . The 
ionospheric delay (X X I is multiple of the Slant Total 

Electron Content (STEC) / , the integrated free electron 
density along the ray that is usually measured in TECU 
units. The term p* represents the terms that do not 
depend on the frequency (range, clock errors, tropospheric 
refraction, etc.). 
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Table 1: The main numbers listed for the different kind of 
observations considered in this work, associated wilh the LI , L2, 
L3, and the extra-wide lane and wide lane carrier phases (Lew 
and Lw): Frequency ( f x . in MHz), wavelength (A Y , in 
meters), Phase Measurement Error (£ Y , in meters), Code 
Measurement Error (E x , in meters), Maximum nnd typical 
Phase Multipath (m x . in meters) , Maximum and typical Code 
Multipath ( M x . in meters), and Ionospheric coefficient {(X x . 
in meters/TEClT). 
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In the case of Wide Area Differential GPS networks, from 
these real-time Slant Total Electron Content (STEC) 
corrections obtained by the tomographic model, it is 
possible, to form the station-satellite double differences, 
VAS7EC , with an error below 0.26TECU, and to obtain 
a second ambiguity (the wide lane) in the reference 
stations. And, finally, to interpolate to the location of the 
rover the unambiguous V&STEC with a precision of a 
few tenths of TECU . If the interpolated value is better 
than 0.26 TECU, then the rover will be able to solve both 
ambiguities in real-time. More details on this technique 
(WARTK) can be found in Colombo et al. (1999) and 
I lernandez-Pajares et al. (2000 a). 

FROM WARTK TO WARTK-3 

We have extend'*.! WARTK, adding the ionospheric 
modelling and correction to the instantaneous three-carrier 
phase ambiguity resolution algorithm (TCAR), to 
medium and long distances (from tens to hundreds of 
kilometers), and with a minimum of additional 
computation. The result is WARTK-3. The main 
improvement is the use of the ionospheric correction in the 
3 rd TCAR step (in previous section). An improvement 
made earlier on to the I st TCAR step will be also 
implemented in WARTK-3: using the mean of the 3 
pseudo-ranges (PI, P2, P3) instead of a single value of 
Pew, to diminish the effect of the multipath in the extra- 
wide lane ambiguity resolution (this would be equivalent, 
in this step, to the ITCAR in single-epoch mode). 
Ionospheric refraction in PI, P2, or P3 has the same sign, 
and nearly the same size as in I^ew, compared widi its long 
wavelength, so this refraction should virtually cancel out 
when subtracting a pseudo-range from the extra-wide lane. 
Finally, an integrity test will be performed using the wide 
lane and LI codes lo delect jumps in the ambiguity 
estimation associated with errors in the extra- wide lane 
ambiguity and wide lane resolution. Indeed, looking at the 
2 nd and 3 rd TCAR steps, equations 5 and 6, it is easy to see 
that 1 cycle of extra-wide lane ambiguity error produces 
an error of about 8 cycles in the regular wide lane and 
each cycle of wide lane error produces about 4 cycles of 
I J error. Many times these ambiguity estimation jumps, 
produced by bad extra- wide! ane ambiguity determinations 
due to the pseudorangc multipath, are large enough to be 
detected and (literal out by using the corresponding Pw 
and PI pseudo-range combinations With three frequencies, 
two additional checks should be made to verify: (a) The 
equal parity (even/odd) of N1-N3 and N1+N3, and of Nl- 
N2 and N1+N2; (b) dial before rounding off lo the nearest 
integer number of wide lane wavelengths, Lew-Pew < 0. 1 
cycle from that integer, and a similar condition on Lw, but 
substituting the pseudo-range with Lc plus its estimated 
bias. These are all "geometry- free" checks on the resolved 
ambiguities. Additional "geometry-dependent" checks 
have been implemented in the long-range kinematic 
software (developed by the last author) that we have used 
in the past, (e.g., Colombo et al., 1999, 2000, 2002): The 
ambiguities are assimilated by the navigation filler as 
pseudo-observations, so the last check is the null- 



hypothesis test made on all a posteriori data residuals to 
either accept or reject the data. As a final precaution, the 
formaf accuracy of a "floated" Lc bias about to be fixed is 
first taken into account if it is better than a few cm, and 
presumably already quite correct, it does riot need fixing, 
and it is not. (The estimation of "floated'' Lc biases, 
whether explicit or implicit, is an integral part of any 
reliable long-range navigation algorithm based primarily 
on carrier phase data.) These "geometry-dependent" tests 
require at least 5 satellites in common view. With only 4 in 
the simulated data sets, we present here results for the 
"geometry- free" part of the technique alone. 
WARTK-3 could help mitigate significantly a major 
limitation on our ability to navigate with high precision at 
long distances (more than 100 km from the nearest 
reference station). And it complements previous 
approaches (see table 2). 





ADVANTAGES 


DISADVANTAGES 


TCAR 


lx>w computational 
load. 


Seriously limited by 
ionospheric refraction. 


ITCAR 


Improved results by 
integrating TCAR in a 
navigation filter. 


The ionospheric delay 
still limits the 3 rd 
ambiguity fixing. 


WARTK 


Accurate real-time 
ionospheric modelling, 
allows precise 
navigation at hundreds 
of kilometers from the 
nearest reference site . 


In spite of speeding- 
up the navigation 
Kalman filter, a 
significant 

convergence lime is 
still needed (5-15 
minutes). 


WARTK-3 


Uses the cxtra- 
widelane, and an 
accurate real-time iono. 
model to provide 
single-epoch precise 
navigation capabilities, 
and greatly speeding 
up the convergence of 
the Navigation Pi Iter to 
just a few epochs). 





Tnble 2: Main advantages and disadvantages of the four real- 
time ambiguity resolution procedures discussed in this work: 
TCAR, Integrated TCAR, WARTK and WARTK-3 frequencies. 



DATA 

To prove the feasibility of WARTK-3, we have analysed 
several dntasets, with simulated airplanes, surface roving 
users, and fixed sites. Wc have studied in detail the 
important case of a receiver on a moving surface vehicle 
(SUR2) at. a distance of e 129 kilometers from a reference 
site (REF5), under different ionospheric conditions. 

The datasets were generated with a modified GNSS 
satellite signal generator, provided by ESA and produced 
in the context of a previous funded study (Laboratory 
Experiment, 2000). The two Cil\S carriers (al 1575.42 and 
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From the last equation, in which additional minor terms 
such as the carrier-phase wind-up are not explicitly 
written, it is possible to estimate the double-differenced 
extra-wide lane ambiguity VA#^ in one epoch, 
substracting the corresponding code: 

VA# W = j- VA(L W - Pj= VA/C - 

In absence of severe multipath affecting the receivers (less 
than few meters), the error term of such estimation, taking 
also into acount the measurement error, is less than 0.5 
cycles (about 3.5 meters), allowing the ambiguity to be 
fixed to the right integer value instantaneously. 

Step 2. The wide lane combination ambiguity (equation 4) 
is estimated from the unambiguous extra-wide lane carrier 
phase obtained in step I. The difference between them 
consists mostly of the wide lane ambiguity, and the 
differential ionospheric refraction (about 0.06 
cycles/TECU with our work frequencies). The non- 
dispersive terms cancel out (equation 5). The main 
problems here are the measurement error and multipath in 
the carrier phase signals. Although typical values of 
differential ionospheric refraction at mid-latitudes and 
baselines below 100 km are just few TECU, the use of an 
ionospheric correction can significantly increase the 
success percentage at longer distances, and in more 
difficult ionospheric scenarios. 

Once the first long-wavelength ambiguities VAA^, are 
solved by using the corresponding pseudoranges, a second 
'mid-wavelenguY ambiguity such as the wide lane 
combination, 

L w =A w( 0 I -0 2 ) = -^-^- = 

= P + KK + a J + m » + e - + - < 4 > 

can be estimated from VAN rw and from the 
corresponding difference of carrier phases (see 
corresponding constant values in table 1): 



VA/9 = VA(L„ - L„ v + X„N„) = VA/V (V + 
K 

- -L V A(e vw + m ew -mj + f(a w - oc cw )V A/ 



(5) 



In the presence of moderate multipath, the corresponding 
error term of the last expression that also includes the 
observational error is usually lower than 0.3 meters, i.e. 
less than 0.3 cycles. The remaining error term of the 
ambiguity resolution corresponds to the ionospheric 
refraction with a value of 0.0580 cycles/TECU (from 
values in table I). This term can decrease the success rate, 
but it is not a critical term at mid-latitudes and at distances 



below few hundreds of kilometers, where the double 
differenced STEC values, VA/ , are typically lower than 
10 TECU (see for instance the corresponding values of the 
analyzed datasets in figure 5). . 

Step 3. The LI phase ambiguity is derived from the 
difference between LI and the unambiguous wide lane 
obtained previously. In this step, the main problem is the 
corresponding ionospheric differential refraction (about 
1.9 cycles/TECU), which can produce errors of several 
cycles at mid-latitudes (see figure 5 with typical values). 

In the third step, we apply a similar approach as in the 
second one, but using the carrier phase differences 
between the short wavelength and the mid- wavelength, 
instead of the mid and long ones: 

VAtf, =^-va(l 1 -l w +a w /vJ=va^ i + 

(6) 

- -1 VA(£„ + m w - m, ) + y- («i ) VA/ 

In this step, the combination of the carrier phase 
measurement error and a moderate multipath only 
introduces again a typical error below 0.2 cycles (see table 
1 ). However, the critical problem here is the ionospheric 
refraction that could introduce errors greater than .5 cycles 
(-1.9475 cycles/TECU, equation 6 and values of table I) 
also for short baselines (see figure 2). This fundamental 
limitation can be overcome by computing real -lime 
ionospheric corrections better than 0.26 TECU (i.e. 0.5 
cycles*! TECU/ 1. 9475 cycles) to assure the right integer 
ambiguity estimation. Once the two wide lanes and LI 
have been resolved, the resolutions of L2 and L3 
immediately follow. 

THE REAL-TIME IONOSPHERIC MODEL 

To overcome the problem in the third step (i.e. to 
determine the shortest wavelength ambiguities), we can 
apply, as in WARTK , a real-time model to estimate the 
ionospheric differential refraction. This model would be 
similar to the one used for meter-level augmentation 
systems such as EGNOS or WAAS, but for sub-decimeter 
carrier-phase navigation it has to meet higher accuracy 
requirements. This model is computed from the dual- 
frequency carrier phase data from fixed network sites, 
assuming a tomographic description of the sounded 
ionospheric region, which reduces significantly the 
mismodelling of the electron content (Hern an dez-Paj ares 
etal. )999a-b). With this approaeh, and the simultaneous 
determination of tropospheric delays and Lc biases with a 
geodetic procedure, it is possible to estimate the 
differential ionospheric refraction with an error below 0.26 
TECU, at distances of hundreds of kilometers from the 
nearest site, under many different ionospheric conditions. 
This precision implies typical errors below 0.5 cycles in 
the LI instantaneous ambiguity fixing at mid- and long 
distances. 
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1227.60 MHz), and the GJjONASS channel 24 carrier (at 
1615.50 MHz) were the three frequencies used in the 
GNSS simulator, with 4 satellites in view for 20 minutes, 
with the data sampling rate of 1Hz, using all 12 available 
channels of an AGGA validation receiver. These data arc 
adequate for our real-time ambiguity determination 
studies, because we are using a geometry-free method, but 
they arc very limited Gust 4 satellites) for navigation 
purposes. Two basic datasets have been considered in this 
work: P5M0 with maximum signal power and no 
mullipath, and P3M1 with medium power and multipath 
(see Laboratory Experiment 2000, AGGA Validation 2000 
for details). In this last, and more difficult situation, the 
corresponding multipath for code and carrier phases, 
mostly affecting steps I and 2 respectively of TCAR, 
reach to more than 4 meters for the pseudo-range, and to 
more than 5mm tor the carrier phase. 

fn addition to the reference receivers simulated in the 
datasets (see layout in figure 4), three additional stations 
(based on existing ones from the 1GS network) have been 
added at distances of more than 200 kilometers to assist 
the ionospheric modelling, emulating a more realistic 
situation with a broader network of fixed sites. 
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Figure 4: Typical baselines from the roving receiver SUR2 to the 
five different reference receivers considered in this work. 

In the original datasct, the ionospheric values 
corresponded to night conditions, with negligible effects 
on baselines shorter than 60 km. In this work, ionospheric 
refraction values corresponding to the noon of the target 
day (17 March 2000, in the Solar Max peak) have been 
generated using an ionospheric ctimatological model (1RI, 
Bilitza 1990). These values are about 5 times greater than 
the Slant Total Electron Content (STEC) values simulated 
in the original data (based on Klobuchar model at night 
time). Other values corresponding to additional scenarios 
(night time, low latitude, high dynamics) have been also 
generated. In figure 5, it can be seen that the STEC, 
double-differenced between pairs of satellites and 
receivers, are larger than 2.7 cm for baselines of 15 
kilometers and more, impeding the third ambiguity 
resolution VAN, (see equation 6). This was not the case 



with the ionospheric values generated in the original, 
datasets, corresponding to night values. Indeed, the double 
differences of the original datasets STEC only exceeded 
the ambiguity resolution threshold when the more distant 
reference stations, at 129 km, is considered. This is 
compatible with the low instantaneous ambiguity 
resolution success rates reported at such distances by 
Vollalhct al.(2001). 

We are going to see in the next section that the application 
of the real-time tomographic mode to correct ionospheric 
delays may often reduce this effect below the 0.26 TEC, 
ensuring the resolution of the ambiguities of all three 
carrier-phases. 



ddSTEC simulated with IRI at noon conditions for the rover SURZ 



REF5-SUF.2 (129 Ho baseline) 



REP1-SUR2 (I baseline) 
>shold to fcoiye ddNl f+/- O.itf TCCU) 




7200 



7400 



7600 7800 0000 8200 H400 B6tM 

GPS TIME (seconds 17 March 2000] 

Figure 5: Double differenced Slant Total Electron Content 
simulated with IRI noon values, VA/ , for surface roving 
receiver SUR2, referred to five different reference stations, 
REFI-5, at distances of I to 129 km far. 



RESULTS WITH THE 
REFERENCE STATION 



MOST DISTANT 



We are going to show the results for the more difficult 
case of the longest baseline: SUR2-REF5, of about 130km 
(see rover trajectory in figure 6). 

Ionospheric filter 

First, let us look at the real-time ionospheric filter 
performance. The estimated versus the true ionospheric 
delays are shown in figure 7. And, more important for 
navigation purposes, in figures 8 and 9, one can see the 
differences (errors) between the double-differenced STEC 
lor the rover according to the model, VA/ , and the true 
values, VA/ , obtained from the rover's data. 92% of the 
values have errors below the threshold limit of 0.26 
TECU, for all the receiver datasets corresponding to this 
baseline. This is accurate enough to be able to resolve the 
3 ambiguities in the absence of multipath, and not taking 
into account the measurement error. Most of the 8% of 
estimates with errors greater than 0.26TECU come from 
satellite PRN 26, that is observed at a low elevation and 
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towards the South, where the highest ionospheric gradients 
take place (see figure 1). These results could be still 
slightly improved by adding ah additional fifth station to 
the south of the present network and of the rover. 



fWr' — FiMPri ttlle REF5 <approx. 1M ka) 
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Figure 6: The trajectory of the ground rover SUR2 is shown. 
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Figure 7: Absolute Slanl ionospheric correction (STEC) 
estimated in real-time with the tomographic model, compared 
with the true values for the surface roving receiver SUR2. 
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Figure 8: Double differenced STF.C est i mated in real-time with 
the tomographic model, compared with the true values for the 
surface roving receiver SUR2 referred to the most distant 
receiver, REF5 (about 129 km away from the rover). 
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Figure 9: Error in the real-time tomographic determination of the 
double differenced STEC. More than 92% of the ionospheric 
estimates present errors below the threshold for successfully 
resolving the third ambiguity VAN, (+/- 0.26 TEC U). 

Ambiguity resolution success 

Once the real-time ionospheric corrections are computed 
and distributed from the reference network to the users, the 
WARTK-3 algorithm can he applied to estimate and fix 
in one single epoch the full set of 3 ambiguities per 
double-di fference. 

A summary of the main ambiguity resolution results is 
shown in tables 3 and 4, for the dalascts P5-M0 and P3- 
Mi, without and with multipath respectively, where the 
success in the three TCAR steps is indicated in three cases: 
(a) without ionospheric corrections, (b) with the 
corresponding Klobuchar model ionospheric corrections 
broadcast by the GPS system and (c) with the real-time 
tomographic ionospheric corrections. 



PS-MO /SUR2-RKK5 
(= 129km) 


Success 


Success 
VA/V\ V 


Success 
VAA/, 


Without ionospheric 
corrections (TCAR) 


100% 


100% 


0% 


TCAR+Klobuchar 
corrections 


100% 


100% 


33% 


Real -time tomographic 
corrections (WARTK-3) 


100% 


100% 


92% 



Table 3: Success relative to the total number of 3834 trials on 
resolving On-Thc-Fly the extra-wide lane, wide lane and LI 

ambiguities (respectively VAA^ t . , VAN V(I and VAN, ) for 

the roving receiver SUR2 referred to the more distant fixed site 
(REFS, at about 129 km ) r and with the ideal datasel P5-M0 
(maximum reception power and no multipath). Each row 
corresponds to cases with different ionospheric corrections: 
without ionospheric correction (first row), with Klobuchar and 
real -lime tomographic corrections (second and third rows). 

It can be seen that the application of the WARTK-3 
procedure, combining TCAR and the real-time accural e 
ionospheric model, allows to increase dramatically the 
single epoch ambiguity determination success not only in 
the ideal scenario P5-M0 (from 0 to 92%) but also in the 
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more difficult one, P3-MI with multipath (tu 92 % f .with 
about 35% using Klpbuchar instead of the tomographic 
corrections). 



P3-M I / SUR2-REF5 
(«I29km) 


Success 
VAN 


Success 

VAN W 


Success 
VAN, 


Without ionospheric 
corrections (TCAR) 


90% 


95% (86) 


3% (2) 


TCAR+Klobuehar 
corrections 


90% 


96% (87) 


35% (31) 


Real-time tomographic 
corrections (WARTK-3) 


90% 


96% (87) 


92% (79) 



Table 4: Same as table I, but for dataset P3-MI (medium level 
of reception power and multipath) and indicating the success 
only for the observations with the TCAR previous ambiguity 
successfully solved (in parenthesis, the success rate regarding the 
total number of observations). 



From table 4, it can be seen that one important problem 
with the P3-MI datasets, affected by multipath, is the lack 
of integrity in the TCAR ambiguity estimations, VAN< W 
and VAN W with about 10% and 4% of erroneous 

instantaneous determinations, respectively. To increase the 
integrity, i.e. to decrease the chance of considering 
erroneous ambiguities, in this case, in die presence of both 
carrier phase and code multipath, we can use the 
corresponding wide lane and LI pseudo-ranges to try to 
detect potential ambiguity determination errors of the 
previous longer wavelength ambiguity in the TCAR 
approach (extra- wide lane and wide lane ambiguities 
respectively). They are amplified by a factor of about 9 
and 4 wavelengths in the extra-wide lane to wide lane and 
wide lane to LI determination, respectively (see equations 
5 and 6). A summary of the corresponding results is 
indicated in tables 5 and 6, using the pseudo-ranges and 
smoothed pseudo-ranges respectively to filter out such 
large errors. 



P3-M1 / 
SUR2-REF5 
(» 129km) 


Success 
VAN * 


Success 

VAN W * 


Success 
VAN, * 


Availability 


Without 
ionospheric 
corrections 

(TCAR) 


99% 


95% (94) 


0(0)% 


38% 


TCAR+ 
Klobuehar 
corrections 


99% 


97% (96) 


33% (32) 


38% 


Real-time 
tomographic 
corrections 
(WARTK-3) 


96% 


96% (92) 


91% (84) 


90% 



Tnblc 5: Similar study 10 thnt summarized in table 4, but {he 
ambiguity success is here computed after passing the integrity 
test described in the text. The last column indicates the 
availability, i.e. % of the 3834 observations passing the previous 
integrity test with the pseudo-ranges (* percentages of 
ambiguities computed after passing the integrity test). 

It can be seen thai with the WARTK-3 approach, a 
significant increase of integrity (from 79% -table 4- to 



91% -table 6-) is attained with a relatively small decrease 
of 16% of availability (from 100% to 84%) using the 
smoothed code integrity check. If the single epoch codes 
are used instead, the integrity is also improved (84%) and 
with a higher availability (90%). Using the GPS Klobuehar 
broadcast model, instead of the real-time tomographic 
corrections, there is a clear worsening of the results, with 
an almost complete lack of availability. 



1WM 1 / 
SUR2-REF5 
(« 129km) 


Success 


Success 

VAN W * 


Success 

VAN, 


Availability 


Without 
ionospheric 
corrections 

(TCAR) 


100% 


100% 
(100) 


0(0)% 


0.2% 


TCAR+ 
Klobuehar 
corrections 


100% 


100% 
(100) 


37% (37) 


0.2% 


Real-time 
tomographic 

corrections 
(WARTK-3) 


100% 


100% 
(100) 


91% (91) 


84% 



Table 6: Similar study to that summarized in table 5, but using 
s moo tlied pseudo-ranges - in the integrity test (* percentages of 
ambiguities computed after passing the integrity test). 



Single- epoch rover navigation 

The datasets available for this study, in spite of being 
adequate to work out the instantaneous ambiguity 
resolution (because this is a "geometry-free" procedure), 
are very limited to determine the single epoch position. 
Indeed, just 4 satellites in view are available during 20 
minutes of data, in such a way that we had to compute the 
instantaneous navigation solution for the surface rover 
SUR2 using the P5-M0 dataset without multipath, and 
with the highest available signal to noise ratio, in order to 
assure the minimum number of 4 satellites for navigation. 
Moreover, the tropospheric delay, orbit errors, etc., for 
which estimation we would need a 5 th satellite at least, has 
been either substracted or left out. Finally the first. 7200 to 
7500 seconds approximately, has been avoided in the 
positioning due to certain carrier phase break patterns. 
They are potentially due to small receiver measurement 
problems in this period. With another dataset containing a 
realistic set of 6 or more satellites (that we hope will be 
available to us in the near future), these hard limitations to 
positioning will disappear, providing better geometry 
(smaller DOP) and the possibility to detect and de- select 
the satellites ( probably at most one, occasionally) 
presenting ambiguity errors that are revealed while using a 
navigation filter, or a Receiver Autonomous Monitoring 
Algorithm, or RAIM (see for instance Parkinson et al. 
1996). 

The corresponding results with the 4 satellites contained 
in the P5-M0 dataset arc indicated in figures 10. They 
contain from top to bottom: The pre- fit residuals, and the 
Fast, North and Vertical components of the instantaneous 
positioning error (all in meters) after using WARTK-3. 
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And, as a reference, the corresponding plots, but with the 
true ambiguities, are shown in figure 1 1. 

The main features of the VAL C pre-fit residuals (top plot 
of figure 10) are the errors around 10 cm associated with 
one cycle error in LI (about 8%, as in table 3). that 
produce single epoch "jump-like" navigation errors, 
represented also in the second, third and fourth plots of 
figure 10, and amplified by the corresponding Dilution Of 
Precision geometric factor (DOP). These bad LI 
ambiguity estimates typically affects one double difference 
of each 3 available per epoch (the minimum number for 
positioning). In such a way that this error affects 3 times 
the positioning (about 24% of the epochs), with 76% of 
time with 3 -D errors below 5 cm and the 100% below 21 
cm. Although the error distribution is not gaussian, the 
resulting 3-D RMS is 7 cm, and 3, 5 and 2 cm for the East, 
North and Vertical components, respectively. 
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Figure 10: Pre-fit double differenced Lc residual (top plot), and 
instantaneous (single epoch) navigation errors in East, North and 
Vertical (up) components (2 nd to 4 th plots) obtained with 
WARTK-3 for the rover SUR2 trajectory at about 130 km from 
the reference site SUR2. 

On the other hand, the erroneous positioning components 
(East, North or Vertical) depend on the geometric 
information provided by the satellite affected by the LI 
one-cycle error. To summarize these results, RMS errors 
of I, I and 2 em are obtained for the Hast, North and Up 
components respectively, when the ambiguity is well fixed 
(about 92% of the trials and 77 % of the epochs) and 3. 5 
and 2 cm of corresponding RMS errors with the full set of 
epochs.EMBED 

The ionospheric free combination Le pre-fit residuals and 
the single-epoch navigation errors, using the true 
ambiguities, can be seen in figure I I. The amplification of 
the carrier phase noise is evident, in particular, in the 
vertical component reaching to about 5 cm of error. This 
trajectory (with and 3-D RMS of 3cm, and 1, I and 2 cm 
in the North, East and Up components, and 95% levels at 
about 2, 2 and 4 cm respectively) represents the solution 
that could be obtained with a real-time filter quite well, 



instead of the single-epoch solution, after Gxing enough 
ambiguities, from the epoch 50 approximately (figure 10). 
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Figure 11: Same as previous figure 10, but for the true 
ambiguities. 

Latency 

In order to characterize the impact of latencies in the 
ionospheric correction (due for instance to the updating 
design of the model, or to potential problems in 
comunications), wc have considered delays from 1 to 30 
seconds in the ionospheric correction computations within 
the fixed sites network. For each of these delays, the 
success percentage on achieving the accuracy of 0.26 
TECU has been computed: the effect is negligible up to 30 
seconds, with success greater than 90%. 
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Figure 12: Similar plot to the previous one, but considering 
common latencies in the reference and rover data.. 

In figure 12, we have also considered greater latencies for 
both reference network corrections and rover in such a 
way that the ionospheric pierce points remain close. It can 
be seen that the previous result is maintained to about 5 
minutes. After this time the success decreases from 90 to 
85%. And after 10 minutes, the success decreases to 75%. 
Note that the jumps correspond to the epochs in which the 
double differenced STEC variations of each satellite reach 
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to 0.26 TECU (see figure 3). These numbers can worsen in 
scenarios with still higher differential ionospheric delay 
variation. These numbers suggest that at mid-latitude the 
ionospheric latency is not an issue for WARTK-3 (up to 
about 5 minutes of latency). 

OTHER SCENARIOS 

In order to better characterize the WARTK-3 performance, 
we summarize the results in different scenarios from the 
point of view of the main problem, the ionosphere, for 
long baselines. The same set of receivers has been 
maintained, but simulating ionospheric nigh I conditions 
(lower values), low latitude (higher values) and high 
dynamics with a roving aircraft (higher variations). 

Night ionospheric conditions 

The performance of WARTK-3 in night conditions with 
small ionospheric delays is still better (100%) than at 
noon time (92%, last section) as it is expected. If the 
ionospheric corrections are not applied, just 62% of 
success is attained in the LI ambiguity determination. 

Lower latitude ionospheric conditions 

In this new scenario, we have considered greater 
ionospheric values, by adopting those corresponding to the 
South of Europe and North of Africa, at a latitude of about 
35 degrees, instead of Central Europe at about 48 degrees 
(figure 1). The success on providing double difference 
ionospheric corrections more accurate than 0.26 TECU is 
still maintained at a relatively high value of about 74%, 
against 33% using the Klobuehar model or just 3% without 
correcting the ionospheric delay. This performance could 
be improved by using the rover ionospheric data as well, 
as in the experiment summarized in the next sub-section. 

Extreme ionospheric con diet ions, under the equatorial 
anomalies 

We present the first results of real-time ambiguity 
determination in the Tropics, below the northern equatorial 
anomaly maximum of the vertical ionospheric delay (see 
figure I). To do this, we have simulated the ionospheric 
conditions at about 22 deg. in latitude over Africa. The 
absolute STEC values and, the double differenced STEC, 
are greater than in the mid-latitude case by about 50- 
100%. To cope with the single-epoch- ambiguity resolution 
at about 130 km from the nearest, reference site in this 
extreme case, we have incorporated, the rover ionospheric 
data, along with that from the stations, as in other studies 
mentioned earlier on. The success attained providing 
accurate ionospheric corrections, after about 10 minutes of 
filter convergence, is around 80% . This result, compared 
to 0% with Klobuehar or with out ionospheric corrections, 
is promising, suggesting the feasibility of" potential 
improvements in future studies. Notice that these results 
are possible because the rover user runs its own 



ionospheric filter, combining its data with the ionospheric 
corrections gathered from the reference neL 
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Figure 13: Map for the high dynamic scenario with the aircraft 
AIR3 flyi ng at the height of 2000 m. 

The last ionospheric scenario studied in this work 
corresponds to the longest baseline at mid-latitude 
conditions, and with the highest dynamics rover available: 
an airplane (AIR3) flying at the height of 2000 m, 
describing ellipses of about 15-20 km in size, about 140 
km away from the reference site, REF5 (see 
corresponding map in figure 13). This cyclic trajectory 
produces wave-like oscillations in the differential 
ionospheric delay (figure 15), which makes this scenario 
more difficult than just with the ground rover SUR2. The 
errors in the real-time double differenced ionospheric 
corrections are represented in figure 16, which 
corresponds to a 73% success for the LI ambiguity, 
against 41% and 5% with Klobuehar and without 
ionospheric corrections. These results could be still 
improved by adding the rover data as in the case of the 
experiment mentioned above. 
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Figure 14: Single epoch estimated STEC and the corresponding 
true values, in the high dynamics scenario. 
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Figure 15: Double differenced STCC instantaneous 
determination, and the true values, in the high dynamics scenario. 
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Figure 16: Instantaneous double differenced STEC error 
in the high dynamics scenario 

RESULTS AT DIFFERENT DISTANCES 

Finally, we present a summary of the results using 
different reference stations, indicated in figure 17. The 
baselines range from 1 km to about 130 km for dataset 
P3M1, and the generated data include multipath effects. 

Indeed, the relative success solving VA/V, (when 
successful wide lane ambiguity occurs), are represented as 
a function of the distance in figure 17, using the datasels 
with medium power reception level, multipath, and IR1 
values lo simulate the ionospheric delays. It can be seen 
that the success decreases from 100% at 1 km to about 
50% at 15 km, and near 0% at 32 km, if no ionospheric 
corrections are considered. However, the inclusion of real- 
time or Klobuchar model estimates of the ionospheric 
corrections maintain the success above the 95% for 
distances of less than 60 km. For the longest bnseiine (at a 
mean distance of 129km), only the tomographic model 
maintains the success above 90%, in contrast with less 
than 40% provided by the Klobuchar model. 
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Figure 17: Relative success of VAN, in the roving 
surface receiver SUR2, as a function of the distance to the 
reference station, computed from the P3-M1 dataset. 

CONCLUSIONS 

The WARTK-3 technique, combining the Three Carrier 
Phase Ambiguity Resolution procedure (TCAR) with the 
real-time tomographic determination of the ionospheric 
corrections, promises the following important 
improvements in ambiguity resolution and precise 
navigation, at mid-latitudes and long distances (more than 
100 km): 

• From an instantaneous (single epoch) success rate of 
about 60% or less in previous works (such as Vollath et 
al. 1998, with night-time conditions), to about 90% at 
distances of more than 100km, and in difficult 
conditions for the ionospheric modelling (noon in Solar 
Maximum). 

• Maintaining high performance even with latencies of 
up to 5 minutes in the ionospheric corrections. 

• The corresponding navigation, in spite of the limited 
number of 4 satellites available in the datasels can be 
done instantaneously, in single-epoch mode, with RMS 
of 3, 5 and 2 cm in East, North and Vertical 
components (about 6, 1 0 and 4 cm at 95%). 

• The good performance has been also checked in other 
difficult ionospheric cases, such as low latitude and 
high dynamics scenarios. 

Additional improvements can be obtained in future work 
by: 

(1) Working with a more realistic, and higher, number 
of satellites, than the four available in the dalasets 
used in this work, would significantly improve the 
navigation quality. 

(2) Using a long-range navigation filter that estimates 
rover position simultaneously with error in the 
tropospheric refraction corrections, errors in the 
broadcast or other predicted GPS orbits, and 
lloated Lc ambiguities as in our fully 
implemented WARTK technique. 

(3) Combining the ionospheric corrections computed 
in the reference network with the rover data in an 
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ionospheric user Filter, such as in the case of 
tropical ionospheric conditions. 

(4) Adding more reference stations for the ionospheric 
determination, to have a better coverage of the sky. 

(5) In cases where it is possible to wait a minute or 
two before the ambiguities are resolved, the 
chances of successful resolution could be 
increased considerably by using the averages of 
several epochs to reduce the effect of noise and 
multipath, as we have done in some of our past 
work with GPS. 



As the main conclusion, the results of this study suggest 
that using a real-time tomographic model of the 
ionosphere, it should be feasible to have instantaneous 
and full ambiguity resolution, achieving sub-decimeter 
precise navigation, at distances of more than 100 km, with 
future three-frequency systems such as GALILEO and the 
modernized GPS. 
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